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Abstract 

Introduction. The widespread use of piezoelectric materials in various industries stimulates the study of their physical 
characteristics and determines the urgency of such research. In this case, modal analysis makes it possible to determine 
the operating frequency and the coefficient of electromechanical coupling of piezoelectric elements of various devices. 
These indicators are of serious theoretical and applied interest. The study was aimed at the development of numerical 
methods for solving the problem of determining resonance frequencies in a system of elastic bodies. To achieve this goal, 
we needed new approaches to the discretization of the problem based on the finite element method and the execution of 
the software implementation of the selected method in C# on the .net platform. Current solutions were created in the 
context of the ACELAN-COMPOS class library. The known methods of solving the generalized eigenvalue problem 
based on matrix inversion are not applicable to large-dimensional matrices. To overcome this limitation, the presented 
scientific work implemented the logic of constructing mass matrices and created software interfaces for exchanging data 
on eigenvalue problems with pre- and postprocessing modules. 

Materials and Methods. A platform was used to implement numerical methods .net and the C# programming language. 
Validation of the research results was carried out through comparing the values found with solutions obtained in well- 
known SAE packages (computer-aided engineering). The created routines were evaluated in terms of performance and 
applicability for large-scale tasks. Numerical experiments were carried out to validate new algorithms in small- 
dimensional problems that were solved by known methods in MATLAB. Next, the approach was tested on tasks with a 
large number of unknowns and taking into account the parallelization of individual operations. To avoid finding the 
inverse matrix, a modified Lanczos method was programmatically implemented. We examined the formats for storing 
matrices in RAM: triplets, CSR, CSC, Skyline. To solve a system of linear algebraic equations (SLAE), an iterative 
symmetric LQ method adapted to these storage formats was used. 

Results. New calculation modules integrated into the class library of the ACELAN-COMPOS complex were developed. 
Calculations were carried out to determine the applicability of various formats for storing sparse matrices in RAM and 
various methods for implementing operations with sparse matrices. The structure of stiffness matrices constructed for the 
same task, but with different renumbering of nodes of a finite element grid, was graphically visualized. In relation to the 
problem of the theory of electroelasticity, data on the time required to perform basic operations with stiffness matrices in 
various storage formats were summarized and presented in the form of a table. It has been established that the renumbering 
of grid nodes gives a significant increase in performance even without changing the internal structure of the matrix in 
memory. Taking into account the objectives of the study, the advantages and weaknesses of known matrix storage formats 
were named. Thus, CSR was optimal when multiplying a matrix by a vector, SKS was optimal when inverting a matrix. 
In problems with the number of unknowns of the order of 10°, iterative methods for solving a generalized eigenvalue 


problem won in speed. The performance of the software implementation of the Lanczos method was evaluated. The 
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contribution of all operations to the total solution time was measured. It has been found that the operation of solving 
SLAE takes up to 95% of the total time of the algorithm. When solving the SLAE by symmetric LQ method, the greatest 
computational costs were needed to multiply the matrix by a vector. To increase the performance of the algorithm, parallelization 
with shared memory was resorted to. When using eight threads, the performance gain increased by 40-50%. 

Discussion and Conclusion. The software modules obtained as part of the scientific work were implemented in the 
ACELAN-COMPOS package. Their performance for model problems with quasi-regular finite element grids was 
estimated. Taking into account the features of the structures of the stiffness and mass matrices obtained through solving 


the generalized eigenvalue problem for an electroelastic body, the preferred methods for their processing were determined. 


Keywords: piezoelectric materials, finite element method, sparse matrices, generalized eigenvalue problem, Lanczos 
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AnHoTayna 

Beedenue. LWinpoxkoe ucnomb30BaHHe Mbe3OMaTepHaJiOB B pa3IMYHbIX OTPacJIAX CTHMYJIMpyeT U3y4eHHe UX (PU3H4eCKHX 
XapaKTepHCTHK HU OOyCOBINBaeT AKTYaJIbHOCTb TaKHX W3bICKAaHHU. B paccMaTpHBaeMOM Cily4yae MOMAJIbHbIM aHasH3 
MO3BOJACT ONpeeIUTh padouylo yacTOTy HU KOIPPUUMEHT ZICKTPOMeXaHHYeCKON CBA3H MIbe3OIIEMCHTOB Pa3JIMYHBIX 
yCTpolcTB. ITH MHAMKATOPh! MpeACTaBIIALOT Cepbe3HbIi TeopeTHYeCKHH UM IpHkaqHOH uHTepec. Len uccneqoBaHua — 
pa3pa0oTKa YHCJICHHBIX MeTOAOB JIA peleHHA 3aa4H OlpeweeHuA YaCTOT pe30HaHca B CHcCTeMe yupyrux Ter. Ja 
OCTWKeHHA We HYKHBbI HOBbIC NOAXOAbI K WHCKpeTH3allMH 3aadH Ha OCHOBeE MeTOa KOHCUHBIX 3JIEMCHTOB 
BBINOJIHeHHe MpOrpaMMHOM peasH3alHu BbIOpaHHoro MeTOa Ha A3bIKe C# Ha IaTPopMe .net. AKTyaJIbHble pewieHHa 
CO3aHbI B KOHTeKcTe OuOMOTeKM KIaccoB KomiIekca ACELAN-COMPOS. OcnHopaHuple Ha oOpallleHHu MaTpHIl 
W3BECTHbI€ MCTOALI pewieHHA OOOOMIeHHOM 3aqa4H Ha COOCTBeHHbIe 3HaYeHHA HeEIPpHMeCHHMBbI K MaTpHuaM OosbUON 
pa3smepHoctu. Jia mpeoxomenua 3TOrTO OrpaHW4ueHHA B pecTaBueHHOH Hay4Ho padoTe peasM30BaHa sJIOrMKa 
MOCTpOCHHA MaTPHL Macc HM CO3aHbl IporpaMMHble MHTepdelich Wid OOMeHa TaHHbIMM O 3aa4ax Ha COOCTBeHHBIe 
3HadeHHA C MOJYIAMH Ipe- MW MOCTIpoLeccuura. 

Mamepuansi u memooo. J\na peanv3alyun 4HCIeHHbIX MeTOROB 3afelicTBOBaIH TaTpopMy .net HU A3bIK 
liporpamMMuposanusa C#. Barnyalua pe3syIbTaTOB HCCIeAOBaHUA NPOBOAMIACh IYTeM CpaBHeHHsA Hal eHHbIX 3Ha4deHHi 
C pelieHwiMu, ouyYeHHbIMH B W3BeCTHBIX CAE-nmakeTax (aHrs. computer-aided engineering — 
KOMIbIOTepH3HpOBaHHasd HHKeHepuA). Co3laHHble MOAMpOrpaMMbl OLCHMBAJIMCh C TOUKH 3peHHA MpOW3BOAMTeIBHOCTU 
MW MIpHMeHHMOCTH Ayia 3aqa4 OombuIoH pasMepHoctH. II[poBoqumMcbh 4YMCIIeHHbIe SKCIEPHMEHTH! C WesbIO Bas Waluu 
HOBBIX aJIFOPHTMOB B 3ajauax MasIOM pa3MePHOCTH, KOTOPbIe pellalOTCA W3BECTHBIMM MeTOAaMu B MATLAB. Jlanee 
TMOAXOA TecCTHpoBaIM Ha 3aqauax C OOJbIIMM UMCIIOM HeM3BECTHBIX HM C yYeTOM pactiapasWiesIMBaHHA OTCJIbHbIX 
onepauui. UroOn u30ex%KaTb Haxo7KeHHA OOpaTHOM MaTPHIbIl, IPOrpaMMHO peasM30BaIH MOAMPUUMpOBaHHbIM MeTOA 


Jlanyoma. Paccmotpemu dbopMaTbl xpaHeHHuA MaTpH B OMepaTHBHOM aMaATH: TpuIIeTE, CSR, CSC, SKyline. [na 
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pelieHHA CHCTeMbI JIMHeMHEIX areOpanyeckux ypaBHenuit (CJIAY) 3ayeticTBoBanM UTepallMOHHbI CHMMeTpHUHbI 
Metog LQ, afanTupoBaHHblii K ITHM (pOpMaTaM xpaHeHHaA. 

Pezyismamot ucciedoesanua. PaspaboTaHbl HOBbIe pacueTHbIe MOY, MHTerpupoBaHHble B OUOMMOTeKy KIaCcoB 
komiiexca ACELAN-COMPOS. IIpopegenbi pacueTsl AIA ompeyeeHHaA IpHMeHHMOCTH pa3IM4HBIX CopMaToB 
XpaHeHHA pa3spe2xKeHHBIX MaTPHI B ONepaTHBHOH MaMATH HU pa3sIM4HbIX MCTOAOB peaslu3alHu Onepalluit c paspexKeHHbIMU 
Matpuiamu. I padwyecku Bu3yaIM3HpoBaHa CTpyKTypa MaTpHL KECTKOCTH, MOCTPOeHHBIX [JIA OHO U TOM xe 3aa4H, 
HO C pa3yIM4HOH TepeHyMepalveli y3I0B KOHe4HOIeMeHTHOM ceTKH. IIpuMeHuTenbHO K 3aqa4ue TeopHu 
3JIEKTPOYIpyrocTH OOOOMIeHEI HU MpeACTaBIeHbI B Bue TAOHMIUbl TaHHbIe O BPeMeHH, HeEOOXOAMMOM Ha BBIMOJIHeHHe 
Oa30BbIX OMepalHii C MaTPHWaMH %KeCTKOCTH B pa3JIM4HbIX (POpMaTax xpaHeHHa. YcTaHOBIeHO, YTO TepeHyMepalHA 
y3J10B CETKH JaeT CyICCTBEHHbIM IpMpOcT MpOu3BOAMTeIbHOCTH JaxKe Oe3 H3MCHEHHA BHYTpeHHeli CTpyKTYpbl MaTpHUbI 
B TaMatTH. C yueTOM MOCTaBJICHHBIX 3aJad MCCIIeqOBaHHA Ha3BaHbl TIPeHMyIMecTBa UM cylaOble CTOPOHbI H3BECTHbIX 
cbopMaToB xpaHeHua MaTpuu. Tax, CSR onTuMasieH Ip yMHOXKeHHM MaTpHIbI Ha BeKTop, SKS — upu obpamenuu 
matpuusl. B 3aqjayax c YMCJIOM HEH3BECTHEIX TWopsKa 10° BEIMTpbIBaIOT B CKOPOCTH UTepaltMOHHble MeTOJbI pelleHus 
oOobmjeHHOH 3aqaun Ha COOCTBeHHbIe 3Ha4eHHA. OUCHHBaIaCh NPOM3BOAUTCIBHOCTb IporpaMMHOH peasM3alHu MeToa 
JIannoma. U3mepsiica Bkay BCex OMepalHii B OOMee BpeMa pelieHHa. BrlicHvsIOcb, YTO ONepalua pemseHua CJIAY 
3aHuMaeT 0 95 % oT oOmjero BpeMeHH paboTsI anropuTMma. IIpu pemenun CJIAY cummMetpwunbim MeToyom LQ 
HavMOOubUIMe BbINMCJIMTeIbHBIC 3aTpaTbl HYKHbI JIA YMHOXKeHHA MaTPHUbI Ha BekTop. Jaq yBemM4eHuA 
MpOH3BOAHTeIbHOCTH aIrOpHTMa UpHOermM K pacnapaswieMBaHHIo c OOmei MamaTBIO. [pH ucnoub30BaHHH BOCbMH 
MOTOKOB IIPOH3BOAUTEIbHOCTh BbIpoca Ha 40-50 %. 

O6cystcoenue u 3akmo4uenue. Ilonyuenuble B paMKax Hay4HOH paooTE! NpOrpaMMHble MOLYJIN OLIN BHEPeHbI B WaKkeT 
ACELAN-COMPOS. OneHeHa UX MTpOH3BOQMTeIbHOCTh JIA MOJCIbHBIX 3aad C KBa3HperyApHbIMU 
KOHCYHORJIEMCHTHBIMH CeTKaMH. C yyeTOM OcOOeHHOCTelM CTpyKTyp MaTPHIl 2%X€CTKOCTH HM Macc, MOyYaeMbIx pu 
pelleHuv OOOOMeHHOH 3aqaun Ha COOCTBEHHBIe 3HAYeHHA JIA 9IEKTPOYMpyroro Tea, OMpeeeHbI pe ANOUTUTeJIbHbIC 


MeTOJIBI JIA UX OOpadoTKH. 


KsroueBble CIOBa: Tbe30MaTepHalibl, MCTO, KOHCYHBIX IICMCHTOB, PaspexKeHHbIe MaTPHUbI, OOOOMIeHHaA 3ayayua Ha 
coOcTBeHHble 3HayeHHA, MeTog JlaHWoma, mosmpoctpanctBo Kppsiopa, MOfYIb lpempoleccuuHra, MOjyJIb 
MOCTIpOWeCCHHTa, TPHIWVIeTbl, KOOPAMHAaTHBIM (bopMaT xpaHeHHaA, CxKAaTHI pa3pexKeHHbIM pad, CSR, cxaroiii 


pa3pexeHHEI crombeu, CSC, SKyline 
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Introduction. Devices made of piezoelectric materials have been widely used, actively studied and improved for a 
long time. Medical ultrasound devices (diagnostic equipment, ultrasonic scalpels) [1-4] and mobile energy generators [5] 
should be noted separately. Paper [6] described the combination of photo- and piezoelectric effects to create efficient 
compact energy sources. New materials designed for the application under specific conditions are being studied in science 
and industry. In [7], the creation of a lead-free piezo-active composition suitable for operation at various temperatures 
was considered. 

In the study on piezoelectric elements, a key role is played by the modal analysis stage, which enables to establish 
the resonance and antiresonance frequencies of the device. These data: 

— are needed to find out the operating frequency of the device; 

— provide determining the electromechanical coupling factor — an important performance indicator of the device; 

— are input information in numerical experiments for problems on forced oscillations. 

The study was aimed at the creation of numerical methods for solving the problem of determining resonance 


frequencies for a system of elastic bodies. Achieving the stated goal required solving two tasks. The first was to develop 
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methods of discretization of the problem based on the finite element method (FEM). The second was to carry out a 
software implementation of the selected method in C# on the -net platform. All known programs take into account the 
context of the ACELAN-COMPOS class library [8]. When solving a generalized eigenvalue problem, methods based 
on matrix inversion are widely used. However, they are not applicable to large-dimensional matrices. In the presented 
research, this limitation was overcome as follows: 

— the logic of constructing mass matrices was additionally implemented; 

— software interfaces were created for exchanging data on eigenvalue tasks with pre- and postprocessing modules. 

Materials and Methods. Principally, the proposed approach was designed to solve static problems of electroelasticity 
when implementing the averaging method [9], which was used to calculate the effective properties of piezo composites. 
In this regard, only stiffness matrices were presented at the stage of constructing global FEM matrices. In this study, we 
additionally implemented the logic of constructing mass matrices and developed software interfaces (application 
programming interface, API) for exchanging data on eigenvalue tasks with pre- and postprocessing modules. The 
developed routines were evaluated in terms of performance and applicability for large-scale tasks. Numerical experiments 
were carried out to validate the algorithms created for such small-dimensional problems that provide obtaining a solution 
by general methods in the MATLAB computing package. Next, testing was performed on tasks with a large number of 
unknowns and taking into account the parallelization of individual operations. 


The mathematical model of the problem being solved consists of the defining relations [9]: 


PreO2U + pPpiau-V-o=$, V-D=0, (1) 

o=cf -(€+ Bs é)-ef -E, D+quD =e;- (€+Ga8) +95 -E, (2) 

6=(Vu+VuT)/2,E=—-Vo. (3) 
Here, o — stress tensor; p; — body density; ¢ — strain tensor; u — displacement vector; D — electric displacement 
vector; E — electric-field vector; f, — body force vector; » — electric potential; a,, By;, Sa — damping 


coefficients; cf, ef 


£. e™, 9% — tensors of elastic constants, piezoelectric modules and dielectric permittivity; index j — 


body number in the model. 
The discretization is performed by replacing: 
u(x,t)=Ni(x)-U(t), ot) = Nj (x): Od). 
Here, N,, — shape function matrix for the displacement field; Ny — shape function vector for electric potential; 


U (t), (rt) — global vectors of the corresponding nodal degrees of freedom. 


In this case, the original problem (1-3) takes the form: 


M-a@+K-a=F. (4) 
Here, matrices M and K are global matrices of mass and stiffness, respectively, and the vector is a general vector of unknowns: 
a=[U, ®]. 
In the problem of the theory of electroelasticity: 
0 K K 
M = M wa ). K _ uu ud . 5 
| 0 0 é uw —K ‘) ” 


Matrices M ws, Kw. and Ky — symmetric. In the case of harmonic oscillations at natural frequency @,, it is 
possible to write: 
a=v,sin(@,t), 
denoting the corresponding eigenvector by 1,. 
Consider free oscillations if F =0. In this case, task (4) is represented as: 
-—07 My, + K-v; =0. (6) 
Thus, the original problem is reduced to a generalized eigenvalue problem (6). For nonzero v,, inequality (6) is solved 


by finding the matrix inverse to K . However, at the same time, the sparse matrix becomes full, i.e., the method is 


unsuitable for large matrices. Therefore, it is needed to use other methods that do not require finding the inverse matrix. 
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To solve this problem, a modified Lanczos method was programmatically implemented in this paper [10]. The author of 
this modification is T. S. Martynova. The description of the development is not given in this article. Of the operations 
used in the method, the most expensive from the point of view of computational resources was the solution to a system 
of linear algebraic equations (SLAE), needed for performing a spectral transformation. 

Matrices M and K — sparse, with a small number of nonzero elements. Several formats are used to store such 
matrices in RAM: 

— triplets, or coordinate format; 

— CSR (compressed sparse row); 

— CSC (compressed sparse column); 

— Skyline storage format (SKS method). 

The coordinate format involves storing triples (triplets) of values (i, j, k), representing coordinates (i, j) and values (k) 
of nonzero elements. CSR is sometimes referred to as CRS or Yale format. It involves storing a sparse matrix in the form 
of three arrays. Consider matrix of N size with NZ nonzero elements. We describe the possible organization of its 
storage. All nonzero elements must be placed in one array of NZ size. The positions of these elements in the columns 
should be placed in another array of NZ size, and the third array of N size should be used to store the indices of the first 
elements of the rows. Similarly, the storage in CSV format is implemented. 

The SKS format assumes the storage of a variable-width matrix band that includes all nonzero elements. In this case, 
zeros are allowed. The efficiency of this format depends on the renumbering of the matrix rows. Methods for reducing 
the size of the tape are described in [11]; however, their applicability to the stiffness matrix obtained when solving a three- 
dimensional problem using FEM requires a separate study. 

To solve the SLAE (system of linear algebraic equations), an iterative symmetric LQ method (SYMLQ [12]) adapted 
to the storage formats listed above was used. 

Research Results. At the beginning of the study, we chose the optimal storage format for sparse matrices. The 
coordinate format enabled to quickly add and change an element of the matrix. These operations were needed at the stage 


of assembling the global matrix and taking into account the boundary conditions. In addition, for ill-conditioned matrices, 


to which K refers, a preliminary transformation is often used for normalization. It is also convenient to perform it in the 
coordinate format. However, this format is ineffective when it comes to algebraic operations. 

CSR is ill-adapted for changing the structure of the matrix: by adding a nonzero element, you need to insert into two 
arrays. In this case, the matrix is multiplied by a vector quite easily and efficiently. 

SKS has similar problems with the addition of nonzero elements and is highly dependent on the renumbering of 
unknowns in the problem. We focus on the example of a quasi-regular grid, which is used in the ACELAN-COMPOS 
package to work with representative volumes of composites. The width of the band containing all nonzero elements can 
be predetermined and depends on the number of nodes and the type of final element. In the general case of an arbitrary 
finite element grid, it is difficult to estimate the size of the band in advance. 

Four methods of numbering unknowns were used in numerical experiments. Figure | shows the structure of stiffness 


matrices constructed for the same task, but with different renumbering of nodes of a finite element grid. 
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Fig. 1. Stiffness matrix structure with various methods of node numbering: 
a — unknowns are ordered by nodes; b — first, displacement nodes are ordered, 
then potentials; c — nodes are sorted by layers of the FE grid, and unknowns — by nodes; 
b, d— nodes are sorted by layers of the FE grid, and unknowns — as in example 


Thus, the grid was a cube with regular partitioning by eight-node finite elements. A model matrix of 500 lines was 
used for the illustrations. The matrices shown in Figures | a and 1 b, were not subjected to additional renumbering of 
nodes and differed only in the numbering of degrees of freedom. In | a: 

a= [Hy My My Py yes Uyy sUyy Uy,» Py | < 

In 1 b, the unknowns responsible for the potential distribution were collected at the end of the vector: 

a= [Uy Myy Uy. yeep Uiyy My Uys > Pyro Py | : 

The unknowns in matrices 1 c and | d were numbered similarly, but the nodes of the finite element grid were pre- 
numbered according to their coordinates through alternately sorting all nodes by each of the coordinates. This technique 
is widely used to build more efficient SLAE solution modules, as it enables to work with the matrix in a suitable band 
format, convenient for parallelization. Similar external modules were implemented for the ACELAN-COMPOS 
complex [13-15]; however, in this work, only formats for storing sparse matrices of a general form were used. 


Table | summarizes the data on the time required to perform basic operations with matrices in various formats. 


Table | 
Time to perform basic operations with the stiffness matrix in the problem of the theory of electroelasticity. 
19,652 rows 
Storage . Elapsed time, ms 
Operation 
format la 1b lec ld 
CSR Conversion from coordinate format 123 132 97 117 
CSR Multiplication by vector, 100 operations 260 260 260 260 
SKS Conversion from coordinate format 690 703 124 268 
SKS Multiplication by vector, 100 operations 60,558 | 61,450 7,616 22,113 


The experimental results showed that the conversion operation from a coordinate storage format to a compact one 
took little time. At the same time, the renumbering of grid nodes to form a block-tape matrix made it possible to get a 
noticeable increase in performance even without changing the internal structure of the matrix in memory. CSR format 
turned out to be optimal in terms of the efficiency of the matrix-vector multiplication operation. When the matrix was 
inverted, SKS format was more efficient, but for problems with the number of unknowns of the order of 103, iterative 
methods for solving a generalized eigenvalue problem worked noticeably faster. 

Further, the performance of the software implementation of the Lanczos method was experimentally evaluated. The 
contribution of all operations to the total solution time was measured. As a result, it was found that the operation of 
solving SLAE took up to 95% of the total time of the algorithm. In the course of the algorithm, a Krylov subspace was 
constructed, and depending on its dimension, the number of SLAE that needed to be solved changed. Note that the 
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dimension of the Krylov subspace was chosen based on heuristics with respect to the number of desired eigenvalues. 
Here, the SLAE differed only in the right part, so that the requirements for allocated memory remained low. Among the 
basic operations used in the course of solving SLAE by the SYMLQ method, the greatest computational costs were needed 
for multiplying the matrix by the vector. 

To increase the performance of the algorithm, the simplest parallelization with shared memory was implemented. 
Blocks of rows were allocated for the CRS format. They were transmitted to separate threads that calculated the 
corresponding components of the resulting vector. The performance gain was 40-50% when using 8 threads. At the same 
time, for matrices of the order of 10° elements, the increase was about 40%, and for matrices of the order of 10*— about 50%. 

Discussion and Conclusions. Within the framework of this study, a method for solving a generalized eigenvalue 
problem for matrices obtained by modeling electroelastic bodies was implemented. Software modules were created in C# 
for constructing mass matrices by the finite element method and performing auxiliary operations within the framework 
of the Lanczos method (working with Krylov subspace vectors, reorthogonalization, finding eigenvectors). The 
computational complexity was mainly due to the operations of multiplying sparse matrices by a vector. In this regard, 
numerical experiments were carried out to determine the optimal formats for storing matrices, the optimal structure of the 
matrix obtained as a result of renumbering the nodes of the FE grid and degrees of freedom in the nodes. A version of the 
SYMMLQ iterative algorithm using parallel computing was developed. The final scheme of work included three points. 
First, global matrices were constructed in coordinate format with a renumbering algorithm (Fig. 1 c). Secondly, the data 
was converted to CRS format. Thirdly, the data was processed by the Lanczos method, which included the SYMLQ 
method for solving SLAE. The results of the work were included in the ACELAN-COMPOS software package. 
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O6 aemopax: 
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MoyemMpoBpaHua MuctutyTa MaTeMaTHKH, MeXaHHKH VM KOMIIbIOTepHBIX HayK «IO2KHbIM PeepabHbIM] YHUBepCuTeT» 
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MaTeMaTHKH, MeXaHHKH MM KOMIIbIOTepHBIX HayK <«FOoKHbI (cbezepalbHbIii yHuBepcuteT» (344006, Pd, 


r. Pocros-Ha-Jlony, ya. bompuraa Cayosaa, 105/42), ORCID, shteyn @sfedu.ru 


Oganesyan PA, et al. Implementation of Basic Operations for Sparse Matrices when Solving a Generalized Eigenvalue 


3aABNEHHbIU 6KIAO Coaemopos: 
ILA. OraHecsH — dbopMupoBaHve OCHOBHOM KOHICHIMH, WesIM HW 3aqja4un UccueqOBaHHA, IporpaMMHad pealm3alua, 
MOTOTOBKa TeKcTa, POPMYJIMPOBAaHHe BEIBOOB. 


0.0. Ute — nporpamMuas peau3alna, IpoBeseHHe YHCJICHHBIX 9KCIeCPHMeHTOB, MOATOTOBKa HIWIOCTpalui. 


Koudauxm uHMepeCcos6. ABTOPHI 3aABIIAIOT 060 OTCYTCTBHU KOH@JIMKTa MHTepecos. 


Bce aéemopbl npOYUmMAaIU U odo6punu OKOHYAMENbHbIU 6apuaHm pykKonucu. 
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